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A DIFFUSION APPROXIMATION THEOREM FOR A NONLINEAR 
PDE WITH APPLICATION TO RANDOM BIREFRINGENT 
OPTICAL FIBERS 1 

By A. de Bouard and M. Gazeau 

CMAP, CNRS and Ecole Polytechnique 

In this article we propose a generalization of the theory of diffu- 
sion approximation for random ODE to a nonlinear system of ran- 
dom Schrodinger equations. This system arises in the study of pulse 
propagation in randomly birefringent optical fibers. We first show 
existence and uniqueness of solutions for the random PDE and the 
limiting equation. We follow the work of Gamier and Marty [Wfroe 
Motion 43 (2006) 544-560], Marty [Problemes devolution en milieux 
aleatoires: Theoremes limites, schemas numeriques et applications en 
optique (2005) Univ. Paul Sabatier], where a linear electric field is 
considered, and we get an asymptotic dynamic for the nonlinear elec- 
tric field. 

1. Introduction. The Manakov PMD equation has been introduced by 
Wai and Menyuk in [31] to study light propagation over long distance in ran- 
dom birefringent optical fibers. Due to the various length scales present in 
this problem, a small parameter e appears in the rescaled equation. Our aim 
in this paper is to prove a diffusion limit theorem for this equation for which 
we will have to generalize the perturbed test function method [5, 20, 24] 
to the case of infinite dimension. In [18, 22], a limit theorem is proved for 
the linear part of the Manakov PMD equation using the Fourier transform 
and the theory of diffusion approximation for random ODE. Obviously the 
method in [18, 22] does not work for a nonlinear PDE. In [12, 22], a limit 
theorem is proved for a nonlinear scalar PDE driven by a one-dimensional 
noise. The proof relies on the fact that the solution processes are continuous 
functions of the noise. These methods are no longer applicable to the limit 
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equation that we will consider which is driven by a three-dimensional noise, 
because the solution cannot be written as a continuous function of the noise. 
Indeed, in a general setting a strong solution of a stochastic equation is only 
a measurable function of the initial data and the Brownian motion driving 
the equation. However, in the case of a one-dimensional noise, Doss [14] and 
Sussman [27] proved that the solution of such an equation can be written 
as a continuous function of the Brownian motion. This result has been ex- 
tended by Yamato [32] to multidimensional Brownian motions when the Lie 
algebra generated by the vector fields of the equation is nilpopent of step p. 
He actually proves the equivalence between the nilpotent hypothesis and 
the fact that the solution can be written as a continuous function of iterated 
Stratonovich integrals. In our case the vector fields driving the Manakov 
PMD equation are functions of the Pauli matrices and the nilpotent hy- 
pothesis of Yamato is not satisfied. This motivates the use of the perturbed 
test function method. Note that the method has been used for a linear PDE 
in [13] and a PDE with bounded diffusion coefficients in [25]. 

We are also interested in the mathematical analysis of both the Manakov 
PMD and the limit equations. Using a unitary transformation, we are able 
to establish Strichartz estimates for the transformed equation, that are not 
available for the Manakov PMD equation. This result will then enable us 
to prove global existence of solutions. The limiting equation is also studied. 
We use a compactness method to study the existence and uniqueness of 
solutions of this latter equation, due to the lack of nilpotent hypothesis and 
to the absence of unitary transformation similar to the Manakov PMD case. 

1.1. Presentation of the model. Optical fibers are thin, transparent and 
flexible fibers along which the light propagates to transmit information over 
long distances and so are of huge interest in modern communications. In a 
perfect fiber, the two transverse components of the electric field are degen- 
erate in the sense that they propagate with the same characteristics: group 
velocity, chromatic dispersion, refractive indices (m = 11,2), etc. However, 
during the fabrication process the fiber may present defects like an elliptic- 
ity of the core or suffer from mechanical distortions like stress constraints or 
twisting [1, 2]. These phenomena induce modal birefringence (n\ ^ 71-2) char- 
acterized by an orientation angle 6 and an amplitude b. If n\ > n-2, we then 
define a slow axe and a rapid axe corresponding, respectively, to the mode 
indices n\ and 77,2- The orientation angle describes the rotation of the local 
polarization axes with respect to the initial axes. The birefringence strength 
(or degree of modal birefringence) is given by b = \m — 7^2 1 ^0 = &i — &2j 
where k\ , ki are the components of the wave vector and ko the wavenum- 
ber of the incident light in vacuum. The beat length Lb = ki-k 2 indicates 
the length required for the polarization to return to its initial state. There 
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exist several types of birefringence that do not have the same effect on the 
electric field. Usually linear birefringence is studied (in the absence of Kerr 
effect, a linearly polarized light remains linearly polarized), although it has 
been shown that the birefringence could also be elliptic (occurring in case 
of twisting, see Menyuk [23]). In case of a uniform anisotropy along the 
fiber, the birefringence parameters (6, b) are constant. However, in realistic 
configurations, the anisotropy is not uniform along the fiber. We assume, as 
in [28-31], that the birefringence is randomly varying, implying polarization 
mode dispersion (PMD). The difference of velocity of the two modes, due to 
random change of the birefringence (and so of the refractive indices) , induces 
coupling between the two polarized modes and pulse spreading: PMD is one 
of the limiting factors of high bit rate transmission. 

In [31], Wai and Menyuk assumed that there is no polarization-dependent 
loss and considered that communication fibers are nearly linearly birefrin- 
gent. We here use one of the models introduced in [31] for which the local 
axes of birefringence are bended with an angle 9 randomly varying along 
the propagation axe and that b and b' (the frequency derivative of b) are 
constant along this axe. Let us recall that the Pauli matrices are defined by 

l\ /0 -A f 1 



1 1 0)> a2= {z 0j> a3 -\0 -1 

and let us consider the coupled nonlinear Schrodinger equation transformed 
into the frame of the local axes of birefringence [21, 31] 



i ■ 



+ £(t)# + ib'a 3 — + 



dt dx 2 dx 2 

+ l\y\ 2 v + i(*v 3 *)a 3 * + In(v) = 0, 

ob 3 

where do is the group velocity dispersion parameter, N(ty) = (^i^l, ^^i)* 
and 

/ b j_<m 

2 dt 

i dO{t) 



V 



-b 



2 dt 

We recall that in the context of fiber optics, x corresponds to the retarded 
time while t corresponds to the distance along the fiber. We introduce a new 
vector field = exp(—ibta^)^ . The evolution of \& is given by the previous 
equation (1.1) replacing S and A r (^), respectively, by 



E(t) 



2 dt 
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Following Wai and Menyuk [21, 29-31], we denote by I the fiber length. We 
also denote by the dispersion length scale and l n \ the nonlinear length 
scale related to Kerr effect. The fiber autocorrelation length l c is the length 
over which two polarization components remain correlated. We consider, 
as in [31], a typical configuration where / ~ Id ~ l ri i ^> l c ^> Lb, that is, we 
consider "relatively small" propagation distances. Under these assumptions 
and the assumptions on d6/dt below, the term N{^) is rapidly oscillating 
and will be neglected [2, 21, 31], its effect being averaged out to zero. As in 
[21, 31], we introduce a unitary matrix 

(1.2) nt)=( M % H\ 

\-U 2 {t) Ml(t) 

the solution of 

(1.3) i^ + E(t)T(t) = 0. 
We also consider, for t S M+, the matrix 

a r u r t)) _ / M 2 " H 2 Vaiu* \_( m 3 m x -im 2 

(1.4) 



2u\u 2 I 1 2 — \ u i\ 2 J \mi + im 2 —7773 
crimi(i) + a 2 m 2 {t) + 037773 (i), 



which characterizes the linear birefringence and where mi, m 2 , 7773 are real- 
valued processes. Then we can remove the rapid variation of the state of po- 
larization in the evolution of ^ using the change of variable *&(t) = T(t)X(t). 
We obtain 

.dX . . ,^.dX d d 2 X 

(1.5) 

+ l\X\ 2 X + -(X*a 3 X)a 3 X + -N U (X) = 0, 

DO D 

where N U (X) = (N 1>U (X), N 2)U (X)Y satisfy 

N liU (X) = (m\ + ml)(2\X 2 \ 2 - |Xi| 2 )Xi 

(1.6) + (772! - im 2 )777 3 (2|Xi| 2 - \X 2 \ 2 )X 2 

+ (mi - irn 2 ) 2 X\X\ + (mi + im 2 )m 3 XfX 2 , 
N 2 , U (X) = (m 2 + m 2 )(2|Xi| 2 - \X 2 \ 2 )X 2 

(1.7) - (mi + im 2 )777 3 (2|X 2 | 2 - [Xil 2 )^ 

- (mi - *m 2 )m 3 Xf Xi + (mi + im 2 ) 2 X 2 X 2 . 

Assuming, as in [18, 22, 30], that the correlation length of dO/dt is much 
shorter than the birefringence beat length and that \d0/dt\ <C b, we set 
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d6 /dt = 2eoa(t), where £o is a small dimensionless parameter and a a Markov 
process with good ergodic properties. Thus, we may replace the process u 
by v, with [18, 22] 

dv(t) = iyJTMivit) o dWi(t) + a 2 v(t) o dW 2 {t)) + i^ s a z v(t) dt 

(1.8) = i^/Tc{<?iv{t) dWi (t) + a 2 v(t) dW 2 (t)) + ij s a 3 u{t) dt 

- i c v(t) dt, 

where |^i(0)| 2 + |^2(0)| 2 = 1, W = (Wi,W 2 ) is a 2d real-valued Brownian 
motion and o denotes the Stratonovich product. The second equation is the 
corresponding Ito equation. In addition, jcjs are two constants determined 
by a and given by 

7c = / cos(26t)E(a(0)a(t)) dt and j s = sin(2&t)E(a(0)a(t)) dt. 
Jo Jo 

Then v(t) € S 3 a.s., the unit sphere in C 2 ~ 1R 4 . We denote by A the unique 
invariant probability measure of v (see Section 5) and by Ea(-) the expec- 
tation with respect to A. Thus, replacing u by v in (1.5), we obtain a new 
equation describing the evolution of the electric field envelope X = (X\, X2) ■ 

6X d &X 8 . 

~ST YThF 9^ ^ 

f)X 1 

= -ib'cr{v(t)) — - -(N V (X) - E A (N U (X))); 

OX 

indeed, the process m = {mi,m 2 ,m 3 ) is now defined as m = (g\ (u), g 2 (i'), 53 (v)) 
and it can be proved (see Section 5) that 

E A (iV v (X)) = |(2|X 2 | 2 - iXxl 2 )^, 
E A (N 2tU (X)) = |(2|^: 1 | 2 - \X 2 \ 2 )X 2 . 

We set 

(1.10) F u{t) {X(t)) = l\X\ 2 X - l(N v (X) -E A (N U (X))). 

Equation (1.9) is of great interest for the study of dispersion because the 
main effects leading to signal distortions (Kerr effect, chromatic dispersion, 
PMD) can be easily identified: on the left-hand side, the first term describes 
the evolution of the pulse along the fiber. The second one corresponds to 
the chromatic dispersion and the last term to the Kerr effect averaged on 
the Poincare sphere. On the right-hand side of the equation, the first term 
describes the linear PMD effect and the second term describes nonlinear 
PMD. 

The Manakov PMD equation (1.9) is written in dimensionless form. Ac- 
cording to the length scales we consider, we set X £ (t,x) = ^A"(^,^) and 
v £ (t) = f(-Ts-), where v is the solution of (1.8); then the electric field X e has 
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the following evolution: 

dx £ (t) ib> dx £ (t) d d 2 x £ (t) (xm _n 

where the term F V£ ^(X £ (t)) is given by (1.10). 

In various physical situations, the long time behavior of a phenomenon 
subject to random perturbations requires to take care of the different char- 
acteristic length scales of the problem. In this context Papanicolaou, Stroock 
and Varadhan [24] and Blankenship and Papanicolaou [5] introduced the ap- 
proximation diffusion theory for random ordinary differential equations. This 
method has been used to study wave propagation in random media [17] and, 
in particular, in randomly birefringent fibers [18, 22], but only few results 
exist on limit theorems for random PDEs. In the latter, the authors studied 
the evolution, in an optical fiber, of the linear field envelope X £ given by 

. ox £ (t) tbf dx £ (t) d d 2 x £ (t) 

and proved that the asymptotic dynamics, when e goes to zero, is given by 

idX(t) + (f ^) dt + ^ o dW k {t) = 0, 

v ' k=l 

where W = (W\, W2, W3) is a 3d Brownian motion, and 7 = (b') 2 /6'j c . Note 
that the linear PMD effect reduces to one single parameter 7 in front of the 
three Brownian motions. Generalizing the perturbed test function method, 
we will prove that the asymptotic dynamic of (1.11) is given by the stochastic 
nonlinear evolution: 

idx { t) + ^^ + F(x(t))y t 

(1.12) +i ^J2°k^-°dW k (t) 

k=l 

= 0, 

where the nonlinear function F reduces to F(X(t)) = ^\X(t)\ 2 X(t) that is 
simply the expectation, with respect to the invariant measure A, of 
F v u\(X e (t)). We will also make use of the following equivalent ltd formula- 
tion: 

((f-f)^+w>)* 

(1.13) +i ^Y^ ak ^l dWk{ t) 

k=l X 

= 0. 
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Note that a different regime concerned with long propagation distances 
and corresponding to I S> l n % ~ 1^ is of physical interest; however, this regime 
would lead to another asymptotic analysis which is beyond the scope of this 
paper. 

This paper is organized as follows: in Section 1.2 we give notation that 
will be used along the paper and state the main results. Section 2 is devoted 
to the proof of well-posedness for the Manakov PMD equation. In Section 3 
we study the local well-posedness of the limiting equation (1.12). Finally, in 
Section 4 we prove the convergence in law of X e to X as e goes to zero. This 
paper ends with Section 5 where we recall some results obtained in [18, 22] 
about the driving process v, and Section 6 where proofs of technical results 
used in Section 4 are gathered. 

1.2. Notation and main results. Before stating the main results of this 
article, let us give some definitions and notation. 

For all p> 1, we define h p (R) = (L P (R;C)) 2 the Lebes gue spaces of func- 
tions with values in C 2 . Identifying C with R 2 , we define a scalar product 
on L 2 (]R) by 



(u, v)u = ^2 J u i y i d%j • 



We denote by W m ' p ,m £ N* ,p £ N* the space of functions in L p such that 
their m first derivatives are in L p . If p = 2, then we denote H m (R) = 
W m ' 2 (M), m £ N. We will also use M~ m the topological dual space of M m and 
denote (•,•) the paring between H m and M~ m . The Fourier transform of a 
tempered distribution v E <S'(R) is either denoted by v or JFv. If s 6 R, then 
M s is the fractional Sobolev space of tempered distributions v G <S'(R) such 
that (1 + \^\ 2 ) s / 2 v £ L 2 . Let (E, \\ ■ \\e) and (F, \\ ■ be two Banach spaces. 
We denote by £(E,F) the space of linear continuous functions from E into 
F, endowed with its natural norm. If / is an interval of R and 1 < p < +oo, 
then LP (I; E) is the space of strongly Lebesgue measurable functions / from 
I into E such that t^ \\f(t)\\ E is in D>(I). The space D>(n,E) is defined 
similarly where (0, J-, P) is a probability space. We denote by L P W (I,E) the 
space L P (I,E) endowed with the weak (or weak star) topology. For a real 
number < a < 1 and p>l, we denote by W a ' p ([0,T], E) the fractional 
Sobolev space of functions u in L p (0,T;E) satisfying 

^ r T \\u(t)-u(s)\\ p 



\ap+l 



— dsdt< +oo. 



'o J o I 1 

The space CP([0,T];E) is the space of Holder continuous functions of or- 
der (5 > with values in E and we denote by Ai (E) the set of probabil- 
ity measures on E, endowed with the topology of the weak convergence 
a(M(E),C b (E)). 
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We will use the space 

x; = (C([o,T],E 1 1 jnc l „([o 1 r],i 1 )nC(o 1 T;H 2 )) xC([o,t],r), 

where C w ([0,T},W n ),m G Z is the space of functions / in L°°(0, T; H m ), 
weakly continuous from [0, T] into HP" - . As the solution of our limit equation 
will not necessary be global in time, we need to introduce a space of explod- 
ing paths, as in [3], by adding a point A, which acts as a cemetery point, at 
infinity in H 1 ; then 

£(M 1 ) = {feC([0,T],M 1 U{A}), 

f(t ) = A for t G [0,T] => f(t) = A for t G [t ,T]}. 

We define a topology on H 1 U {A} such that the open sets of H 1 U {A} are the 
open sets of H 1 and the complementary in H 1 U {A} of the closed bounded 
sets in H 1 . For any / G C([0,T],H 1 U {A}) we denote the blowing-up time 
r(f) by 

r(/)=inf{tG[0,T],/(t) = A} 

with the convention r(f) = +oo if f(t) ^ A for all t G [0, T]. We endow the 
space ^(IHI 1 ) with the topology induced by the uniform convergence in H 1 
on every compact set of [0, r(/)). 

Let (A,G,Q) be a probability space endowed with the complete filtration 
(Qt)t>o generated by a two-dimensional Brownian motion W = 
which is driving the diffusion process v given by (1.8). We first state an 
existence and uniqueness result for (1.11). 

Theorem 1.1. Let e > and suppose that X £ (0) = v£ L 2 (BL), then there 
exists a unique global solution X £ to (1.11) such that, Q^-almost surely, 

X £ G C(M+,L 2 ) n C^M^HT 2 ) nLf oc (M+,L 4 ). 

Moreover, equation (1.11) preserves the L 2 norm, that is, for all t G M+ 

ll^e(*)lk 2 = 1Mb- 

If, in addition, X £ (0) = w£H 1 (resp., H 2 , resp., H 3 ), then the corresponding 
solution is in C(R + ,M 1 ) [resp., C(R+,M 2 ), resp., C(M+, H 3 )] . 

Let (0, J-, P) be a probability space on which is defined a three-dimensional 
real-valued Brownian motion W = (W\,W2, W3). We denote by (J 7 t)teM.+ the 
complete filtration generated by W. The next theorem gives existence and 
uniqueness of the local solution for (1.12) 

Theorem 1.2. Let Xq = v G 1H 1 (1R), then there exists a maximal stop- 
ping time t*(v,oj) and a unique strong solution X (in the probabilistic sense) 
to (1.12), such that X G C([0,r*),H 1 (R)) F-a.s. Furthermore, thel? norm 
is almost surely preserved, that is, Vt G [0, r*), ||A"(t)|| L 2 = \\v\\i2 and the fol- 
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lowing alternative holds for the maximal existence time of the solution: 
T*{v,uj) = +00 or limsup ||X(i)|| H i = +00. 

Moreover, ifveW 2 , then X E C([0, r*), M 2 (R)) and t* satisfies 
(1.14) t*(v,u) = +oo or lim ||X(t)|| H i = +00. 

t/*T*(t>,Ul) 

Note that we do not obtain global existence for (1.12), due to the lack of 
control of the evolution of the H 1 norm (see Remark 3.1). 

Using these existence theorems, we are able to prove a diffusion approxi- 
mation result for the nonlinear system of PDEs (1.11). 

Theorem 1.3. Let X e (0) = X = v be in M 3 (R), then the solution X £ 
of (1.11) given by Theorem 1.1 converges in law to the solution X of (1.12) 
in £ (H 1 ), that is, for all functions f in C&(£(H )), 

lim £(X £ )(/) = £(*)(/). 

Note that we consider here the Manakov PMD equation (1.11), but the 
method may be carried out to other nonlinear Schrodinger equations. Let 
us first emphasize the key points that allow us to prove Theorem 1.3. 

The first point is that the noise term is a linear function of the unknown 
X £ . This particular structure leads to a stochastic partial differential equa- 
tion for the limiting equation. The second point is the fact that the Pauli 
matrices are Hermitian. This is important to obtain the conservation of the 
L 2 norm for both equations. Finally, we use that the driving process v is a 
homogeneous Markov ergodic process defined on a compact state space such 
that Ea(<t(t/)) = 0. The hypothesis on the driving noise may be weakened 
as in the case of a random ordinary differential equation assuming good 
mixing properties (e.g., exponential decay of the covariance function). The 
boundedness of a(v £ (t)) seems to be necessary. It is used to prove uniform 
bounds in Lemma 4.5 for tightness. On the other hand, the lack of Strichartz 
estimates for the limiting equation (1.12) is a negative aspect. Thus, we use 
that F(v) is locally Lipschitz in H 1 (M) to prove existence and uniqueness 
of a local solution to (1.12). But if a{v £ {t)) were a one-dimensional process, 
larger dimension and larger power in the nonlinear term could be considered. 

Other types of nonlinear Schrodinger equations may be considered re- 
placing, for example, i by X £ and assuming that the matrices Of, are 
real valued and symmetric. This latter equation is simpler to handle using 
Strichartz estimates for the fundamental solution and because a{v £ {t))X £ {t) 
can be treated as a perturbation as far as we are concerned with existence 
of solutions. 
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2. The Manakov PMD equation: Proof of Theorem 1.1. The point here 
is that no Strichartz estimates are available for (1.11) because of the lack 
of commutativity of the matrix a at a different time: a{v(t))<r(v(s)) ^ 
a(v(s))(r{v{t)). Consequently, only local existence and uniqueness for initial 
data in H 1 can be easily proved directly on (1.11). The idea of the proof 
is then to find a unitary transformation such that Strichartz estimates are 
available for the transformed equation. This change of unknown is given in 
the next result. 

Lemma 2.1. Let us denote for t G M + 



where v e = v{t/e 2 ), v given by (1-8). Assuming that X e G C([0, T], L 2 ), we 
set ^> £ (t) = Z e (t)X e (t) ; then the evolution of the electric field ^f e is given by 
the stochastic ltd equation 



(2.1) 



. ib' <9# £ d d 2 ^ £ 5 IT|2t 1 /Tllt _ , T 

^ ~di + f ~tW + G 1 ^ 1 ^ + 6 ^e^e)o^ e 

+ ^<j^ £ dt + dt - ^(a^ £ dW x {t) + a 2 ^ £ dW 2 {t)) = 0, 

where Wj(t) = £Wj(t/£ 2 ),j = 1,2, and with initial conditions 

^ (0) -i-, 2 , £ (0), 1+ I7 lie (0),J-^ 

Proof. Using the equation satisfied by v £ and because |^i >£ (£)| 2 + 
I ^2,6 (*) 1 2 = 1 f° r an y ^ — 0' we obtain 

idZM)Z7 x ^ e (t) = ~a 3 ^ £ dt--^^ E dt + ^a 1 ^ £ dWi(t) 

£ z £ z £ 

+ ^a 2 V £ dW 2 (t). 
The nonlinear part of (2.1) is obtained as in the derivation of (1.5). □ 
We first investigate the behavior of the linear equation 

(2 - 2) i ^r + £ iba3 ^ + Y^ =0 

with initial condition ^(O) = tpo G L 2 . 

Proposition 2.1. The unbounded matrix operator H e = ^h-^z — 

e" " 3 ~§x defined on S#{H £ ) = H 2 is the infinitesimal generator of a unique 
strongly continuous unitary group U £ [t) on L 2 . Moreover, U e (t) may be ex- 
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pressed as a convolution kernel, that is, for -ipo G <S(R) 
U £ (t)ip = A £ {t)*^o 



1 



/ [iix-b'tlef 
/eXP \2 dot 



^lixid^t 







V 







ex Pi -- 



i (x + Vt/e) 



dot 



PROOF. Assuming tjjQ G <S(R) and taking the Fourier transform, in the 
space variable, of (2.2), we obtain readily 



1-,/ 



at e" ° " 
Since o~ 3 does not depend on time, we obtain 



* s (t) = R E (t)^ 



( J ^0 

' exp< 



idn 



.6' 



V 



i A t-i-it 

2 S £ 




exp 







^0- 



The statement of Proposition 2.1 follows then in a classical way, setting 
A £ (t)=T~ 1 (R £ (t)). □ 

The explicit formulation of the kernel given in Proposition 2.1 allows 
immediately to get the following dispersive estimates: if p > 2, t 7^ 0, then 
f7 e G £(L P ', L p ) where p' is such that 1 + £ = 1 and for all ^ G l/, 

(2.3) ||^(t)Vo||L P < (2vr|d ||t|)" 1/2+1/p ||V'o|| L p'. 

Using then classical arguments (see [7, 19]), one may prove Strichartz in- 
equalities for U £ (t). 

Proposition 2.2. The following properties hold: 

(1) For every ^ G L 2 (R), U e (-)if) G L 8 (R;L 4 ) n C(R;L 2 ). Furthermore, 
there exists a constant C such that 

I|£4(-)^o||l8(R;L 4 ) < C||^o||l 2 for every ip G L 2 . 

(2) Let I be an interval of R and i G d. det / G L 8 / 7 (d, L 4 / 3 ), then the 
function 



t^ / U £ (t-s)f(s)ds 



to 



belongs to L 8 (d,L 4 ) n C(d,L 2 ). Furthermore, there exists a constant C in- 
dependent of I such that for every /GL 8 / 7 (d,L 4 / 3 ) 



U £ {--s)f(s)ds 



to 



L 8 (/,L 4 )nL°°(/,L 2 ) 



< C|l/llL8/7(/,L 4 /3)- 
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We now turn to the study of the nonlinear problem. We will use, as is 
classical, a cutoff argument on the nonlinear term which is not Lipschitz. The 
cutoff we consider here is of the same form as the one considered in [9] . We 
first prove an existence and uniqueness result for this truncated equation, 
then deduce from this result the existence of a unique solution for (2.1). We 
denote: 

/(* e ) = ||^| 2 ^ + i(**a 3 * £ )^3* £ . 

Let 9 G C£°(R) with suppG C [-2; 2] such that @(x) = 1 for |x| < 1 and 
< ®(x) < 1 for x £R. Let R > and @r(x) = @(x/R). We then consider 
the following equation: 

*f (i) = U B (t)ik + ^ / U E (t-s)a 3 t??(s)ds 
£ Jo 

f U e (t-S)*?(8)d8 

£ Jo 

(2-4) +i [ U £ (t-s)@ R (\\^\\ L s iQ ^ ) )f^f(s))ds 

Jo 

■ I — ft 

_ hHl / U £ (t-s)a l ^(s)dW l (s) 
e Jo 

■ I — ft 

- hLk / U £ (t-s)a 2 y?(s)dW 2 (s), 
e Jo 

which is the mild form of the ltd equation, 

(2.5) -^a 1 ^dW 1 (t)-^a 2 ^dW 2 (t) 

+ e R (\\^\\ L s {Q>t . M) )f^(t))dt = 
with initial condition ^^(0) = f/V 

Proposition 2.3. Let^f (0) = € L 2 (R). LetT>0 andUj = C([0,T]; 
L 2 ) ni 8 (0,T;L 4 ); £/ien f^.^J /ias a unique strong adapted solution £ 
L 8 (A;Uj), for anyT>0. 

Proof. We use a fixed point argument in the Banach space L 8 {A\lAj) 
for sufficiently small time T depending on R. We first need to establish 
estimates on the stochastic integrals 

JjMt)= I U £ (t- 3)^^(8) dWj(s), J = 1,2. 
JO 
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Lemma 2.2. Let T > 0; then for each adapted process ^ £ G L 8 (A;U^) 
and for j = 1,2 the stochastic integral Jj^e belongs to L 8 {A]Uj) . More- 
over, for any T > and t in [0,T] we have the estimates 

^GI^e^elli^O.T^nL^^TjL 2 )) - CT 4 E(|| * e |l!°o( 0>T;]L 2) )• 

Proof. Since ^ £ G L 8 (A;U]:) and is adapted, we may apply the Burk- 
holder-Davis-Gundy inequality in the Banach space L 4 (M) (which is UMD 
space [6]): 



E (ll«^,£*e|lz,8(o,T;L 4 )) 



E 



T 



U £ {t-s)a j ^ £ {s)dW J {s) 



dt 



L 4 



< 



< 



T 



E{ sup 

vO<n<t 



U £ {t-s)a^ £ {s)dW j {s) 



L 4 



dt 



C®(J (J \\U £ (t-s)a^ £ (s)\\l i ds^ dtj. 
Using the Holder inequality in time, Fubini and a change of variable, 

E (i) Oo l|C/e( * _S)CTj '* £(s)l| L 4dS ) dt 



o 



|E/ e (-)oj-* e (s)|||g (0iT . L 4 ) ds 



On the other hand, by Proposition 2.2, 

||C/ £ (-)a- i ^ £ (s)||| 8(0iT . L 4 ) ds 



o 







l= 2 d s 



<CTE(||tf " 8 



; £llL°°(0,T;L 2 )V- 

Combining these inequalities leads to the estimate in L 8 (0,T;L 4 ). The other 
estimate is proved using the Burkholder inequality in Hilbert space and the 
unitary property of the group U £ . Finally, U £ (t) being a unitary semigroup 
in L 2 , Theorem 6.10 in [8] tells us that, provided * e G L 8 (A, L 2 (0, T; L 2 )), 
then Jj £ ^ r e(-) has continuous modification with values in L 2 (M). □ 

Given *f G L 8 (A;Uj), we denote by T^f(t) the right-hand side of (2.4). 
Since the group U £ {-) maps L 2 (M) into C(M,L 2 (M)), Proposition 2.2 and 
Lemma 2.2 easily imply that the mapping T maps L 8 {A;IA]:) into itself. 
Let now and $f being adapted processes with values in L 8 (A;U^), 
then using Proposition 2.2, the same arguments as in [9] for the cutoff and 
Lemma 2.2 applied to Jj i£ ($f{t) - *f (*)), we get 

E(||T*f - T*f ||^t) 1/8 < + + C( j R)T 1 / 2 ')e(||^ - ^|| 8 t) 1/8 . 
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We conclude that T is a contraction mapping if T is chosen such that 
CT/e 2 + CT l / 2 /e + C(R)T 1 / 2 < 1. As usual, iterating the procedure, we de- 
duce the existence of a unique solution of (2.4) in L 8 (A;Uj) for all T > 0. □ 

Our aim is now to get global existence for the process *& e , the solution of 
(2.1) which may be constructed from the above results. Let us set 

= mi{t > 0, 11*^11^8(0,^4) > R}, 

which is a Q £ {t) stopping time. It can be proved using Strichartz estimates 
and the integral formulation (2.4) (see [9, 10]) that is nondecreasing 
with R and that *f = tyf on [0,«f] for R < R' . Thus, we are able to 
define a local solution * e to (2.1) on the random interval [0, where 
K*(ipo) = limR-!.+oo by setting * e (t) = ^f{t) on [0, nf]. It remains to 
prove that k* = +oo almost surely. From the construction of the stopping 
time k* it is clear that a.s., 

(2.6) if K* £ (ip ) < +oo then lira J^fW^fo^iA) = +°°- 

*/*KjW-o) 

The arguments are adapted from [9]. We first prove the following lemma: 

Lemma 2.3. Let ^ e (0) = tpo be as in Proposition 2.3 and be the 
corresponding solution of (2.5); then for any t <T 

\\^s(t)\y = Uo\y a.s., 

and there is a constant M £ > 0, depending on T and ||^o||l 2 ; but independent 
of R, such that 

(2-7) E(II*?IIl8 (0 ,t;L4)) <M £ (T). 

Proof. To prove that the L 2 norm of the solution of (2.5) is constant 
in time, we apply formally the ltd formula to ^||*^(t)|| 2 2 and notice that 
by integration by parts 

tin S *? fa" tin 8 *?"\ n 

Since a* = aj,j = 1,2,3, where * stands for the conjuguate transpose, we 
get 

(vtf (t),i*j*? =0 for j = 1,2,3. 
Moreover, because the ltd corrections cancel with the damping term — jrf^f 
of (2.5), we get ||\E'^(i)|| L 2 = ||-0o||l 2 > ^ ^ T. The computations can be made 
rigorous by a regularization procedure. 

In order to prove (2.7), we follow the procedure in [9, 10]. Using the inte- 
gral formulation (2.4), the conservation of the L 2 -norm and Proposition 2.2, 
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we obtain for a.e. wGl! and for all time T\ such that T > T\ > 

(2.8) ||*f|| L 8 ( o,T 1 ;L*)<^H + CTl /2 ||*fll|8(0,T 1 ;L 4 )» 

where 

K e (u) = Cfl + ^) ||Vo||l2 + - Yl I' J j,e^e\\L^0.T^)- 

v £ J £ j=1 

From inequality (2.8) it follows that H^eIIls^Ti;!/) — 2if e (u)) if Ti is chosen, 
for example, such that Ti(uj) = inf(T, 2~ 6 (C 1 / 2 K £ )~ 4 ). If T x < T we can 
reiterate the process on small time intervals [IT±, (I + l)Ti] C [0,T] (keeping 
R fixed and varying I) to get || II z, 8 (zt i i ;L 4 ) < 2K e (pj). Summing these 
estimates, using Ti = 2- 6 C- 2 (K e y* and the Young inequality, we obtain 

\\*?\\ L *(o,T^)<C(T)(K £ (u;)f. 

Taking the expectation in the above inequality, using the Holder inequality 
and Lemma 2.2, we get the following estimate: 

(2.9) E(||^|| i8(0>T . L4) ) < C(T)((l + ^YwoWl + ^Hl^) , 

from which (2.7) follows. □ 

We easily deduce from Lemma 2.3 and (2.6) that k* = +oo a.s. and as 
in [9] the existence and uniqueness of a solution \& e of (2.1), a.s. in for 
any T > 0. 

To end the proof of Theorem 1.1, we have to extend those results to the 
process X £ . For a.e. u in A and for each t > we set X e (t) = Z~ l (t)^ £ (t) . 
By definition of the process Z~ 1 (t) [which, in particular, is measurable with 
respect to Q £ {t)} and properties of we easily deduce that X e (t) is adapted 
and continuous with values in L 2 , and satisfy (1.11), hence is C 1 with values 
in H -2 . By unitarity of Z £ we also deduce that for all t > 

||* e (t)|| 2 2 = (X e (t),Z-\t)Z E (t)X £ (t)) h 2 = \\x £ (t)\\h, 

and since the coefficients of Z~ 1 (t) are a.s. uniformly bounded, X £ £ L 8 oc (R + , 
L 4 ) a.s.; Theorem 1.1 is proved. 

We now extend the previous global existence results to more regular initial 
data. T being fixed, we denote 

V T = L°°(0, T; M 1 ) n L 8 (0, T; W 1 ' 4 ) 

and 

Vj = C(0, T; H 1 ) n L 8 (0, T; W 1 ' 4 ). 

PROPOSITION 2.4. Let * e (0) = -00 G H 1 and let T > 0; then equation 
(2.1) has a unique strong solution ^ £ with trajectories in C(0,T;M 1 ) . 
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PROOF. Let be in H 1 . Given $f G L S (A; V T ), we denote by (t) 
the right-hand side of (2.4) and W r = L°°(0, T;L 2 ) n L 8 (0, T; L 4 ). By Propo- 
sition 2.2, Lemma 2.2 applied to d x ^^ and the Holder inequality, we deduce 
that 

\\Td x ^\\ L&{uT) < C\\d x Mv + (^T + + CTWuA \\d x ^\\ L s iU Ty 

Therefore, we conclude that choosing Rq = 2C|| v I / o||h 1 i T maps the closed 
ball of L 8 (A;V T ) with radius Rq into itself, provided T is small enough 
depending only on R and e, but not on Rq. Combining with the fact that T 
is a contraction in L S (A;U T ) and that the balls of L 8 (A;V T ) are closed for 
the norm in L 8 (A;U T ), we conclude to the existence of a unique fixed point 
q,R £ L 8 (A; V T ). Using Proposition 2.2 and Lemma 2.2, we get continuity 
of the solution in H 1 . Since the cutoff only depends on the L 8 (0,T, L 4 (]R)) 
norm, we deduce that there is a unique global solution ^ £ to (2.1) with 
paths in C([0, T]; H 1 ). Since the transformation Z £ does not depend on x, 
we conclude that these results still hold true for X e . □ 

Proposition 2.5. Let * e (0) =^ £ M m , m = 2,3. Let T > 0; then 
equation (2.1) has a unique strong solution ^f e with paths in C([0, T]; M m ), 
m = 2,3. 

Proof. We consider equation (2.5) but with ©rGI^^IIls^^l 4 )) replaced 
by 6^(11^^)11^). Given «f in L 8 (A; L°°(0, T; H 2 (M))), we denote by 
T^f^(t) the right-hand side of the integral formulation of this equation. 
We easily prove that T maps the closed ball of L 8 (A; L°°(0, T; H 2 (R))) with 
radius Rq into itself, for Ro = 2C||\E f o||H 2 > provided that T is small enough, 
depending only on R and e, but not on Rq. Using that this ball is closed 
for the norm in L 8 (A; L oo (0,T;IHI 1 (IR))) and that T is a contraction for the 
norm in L 8 (^l;L oo (0,T;IHI 1 (lR))), we deduce that there exists a unique solu- 
tion ^ £ with paths in C(0, T; H 2 (1R)) a.s., which is global since the solution 
is global in Mr. Existence and uniqueness in H 3 can be proved by the same 
arguments. Again those results are easily extended to X e and this concludes 
the proof of Theorem 1.1. □ 

3. The limiting equation: Proof of Theorem 1.2. In order to prove a 
local existence and uniqueness result for the system (1.12), we use a com- 
pactness approach (see, e.g., [16]) motivated by the fact that we do not 
know if Strichartz estimates are available for (1.12). Indeed, no transfor- 
mation similar to the Manakov PMD case seems to be available, as the 
equation dX(t) = — ^/7^fc=i a k dWj~(t) cannot be solved in a simple 
way. We first prove existence of a unique solution in H 1 for the linear part 
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of the equation, defining then a random propagator, and then consider the 
nonlinear part as a perturbation. We will strongly use the fact that the non- 
linearity is locally Lipschitz in H 1 . The regularity in H 2 will follow with the 
same arguments as for (2.1). Let us consider the linear part of (1.12), 

«M = (40)*-^X>^lr°<^M 

(3.1) 

/do 3 7 \3 2 I, 8X(t) mr „ 

v 7 fc=i 

with initial data X(0) = v£ H 2 . We introduce, for 7j > 0, the mollifier = 
{1 — 7) ^-) _1 . We denote by X v the solution of the regularized Ito equation 

(3.2) = (if + 1) ^ * - v=r S> -WiM 

and X^(0) = t)6l 2 . Since the operators <9 2 J 2 and <9 X are bounded from H 1 
into H 1 (with constants depending on 77), we easily get, thanks to the Doob 
inequality, the Fubini theorem, the Ito isometry and the independence of 
(Wfc)fc=i ; 2,3, the existence and uniqueness of a solution X v to (3.2) with paths 
in C([0, T], H 2 ) for any T > 0. Moreover, it is easy to see that the H 2 norm 
of Xjj is conserved since the Pauli matrices are Hermitian. Consequently, 
the process 

f l {idn 37 \ d 2 JlX n 
M,{t) = -X v (t) + X„(0) + J o {-Y + i) ds 

is a JFt martingale with paths in C([0,T],L 2 ). Let us compute the quadratic 
variation. Let a = (ai, 02)* and 6 = (61, 62)* be in L 2 and T > t > s > 0; then 

E((a,M^(t)) L2 (6,M^(i)) L 2 - (^M^^M^s))^) 

We deduce that the quadratic variation of M„(t) is given by 

(3.3) (b,({MS)))a)^=jj2 f( a ^^^) (^k^A du. 

f^Jo V dx J L2 \ dx J L2 

Using the conservation of the H 2 norm and equation (3.2), we get for all 
0<a<± 

(3-4) E(||Xj Ca([0)T];L2) )<C a (T), 
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where C a (T) is a constant independent of r\. Using the Ascoli-Arzela and 
Banach-Alaoglu theorems, the Markov inequality and inequality (3.4), we 
get that the sequence (£(X V )) V>0 is tight on C w ([0,T],M 1 (M))nL2 D (0,T,lH 2 ). 
The Skorokhod theorem [4, 15] implies that on some probability space 
(0, J 7 , J-"t,P), there exist a sequence of stochastic processes (X v ) v> q, and 
a process X, such that 

C(X V ) = C(X V ), C(X) = C(X) 

and lim^o X V = X, P-a.s. in C w ([0, T], M 1 ) n L%(0, T, M 2 ). For all r) > and 
t £ [0,T] we define the process 

m = -x^) + x vi o) + jT (f + f ) *, 

We deduce from the above laws equality that M n {t) is a square integrable 
continuous martingale with values in L 2 with respect to the filtration J~t and 
that the quadratic variation ((M v (t))) is given by formula (3.3) replacing X v 
by X„. Let a G H , then by the above martingale property we get for all 

s < t 

E((a,M v (t)-M r) (s))^\T s )=0. 

Using the almost sure convergence in C lu ([0,r],BI 1 (]R)) of X v , the bounded- 
ness in H _1 of the operator J v and the conservation of the H 1 norm, we get 
the almost sure convergence in C u ,([0,T],]HI _1 (R)) of M„ to M, where 

M (1 )^( t )-i(»)-f(f^)f^. 

Hence, M is a weakly continuous martingale with values in HI -1 . Moreover, 
using the a.s. convergence in C W ([0,T],I1 1 (R)) and the dominated conver- 
gence theorem, we get for all t, s 6 [0, T], t > s and for any a, b E H 1 , 



UmE((6,«M„(t)»a)|^) 



Thus, the quadratic variation (6, ((M(t)))a) is given, for all t G [0, T], by 
(3.5) (b, «M(t)»a> =-vYlJ*( a,(7k ^ {u) )( b,ak l^ {u) ) dU - 
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Noticing that M (0) = and using the representation theorem for continuous 
square integrable martingales, we obtain that, on a possibly enlarged space 

(fi, J-, J-"t,P), one can find a Brownian motion W = (W\, W2, W3) such that 
(a,M(t)) = V7 fY,( a ^^{s))dW k {s). 



Thus, we deduce that (X,W) is a weak solution of (3.1) on (fi, J 7 , .Ft,P) 
with values in C^QO^H^R)) n L°°(0,T,H 2 ). To conclude the proof, we 
have to prove pathwise uniqueness of the solution and strong continuity 
in H 1 . Since X G L°°(0, T, HI 2 ) is the solution of (3.1), we easily deduce 
that X G C a ([0,T],L 2 ) for any a G [0, 1/2). By interpolation we obtain that 
X G C([0, T], HI 1 ). It follows, using the ltd formula, that pathwise unique- 
ness holds for (3.1) in CQOjT],!! 1 ). This implies, by the Yamada-Watanabe 
theorem, that the solution exists in the strong sense. Thus, we can define 
a random unitary propagator U(t, s) which is strongly continuous from H 2 
into H . This random propagator can be extended to a random propagator 
from H 1 into H 1 using the continuity of X in FJ 1 , the density of FJ 2 into H 1 
and the isometry property of U(t,s) in H . 

The local existence of the nonlinear problem (1.12) in H 1 follows from 
the construction of the random propagator U : we consider a cutoff function 
9 G C£°(R), 9 > satisfying 

1, if \\X(t)f ml <R, 
0, if ||X(t)|| 2 1 >2 J R, 

and first construct a solution X R of the cutoff equation, 



®R(\\X(t) 



I Ml , 



idX R (t)+ ^^ + G ii (||X R (t)|| 2 1 )^ ii )wj dt 

(3.6) + E ^ dW k (t) 

k=l 

= 

with initial data X R (0) = »£i' and whose integral formulation is given 
a.e. by 

(3.7) X R (t) = U(t,0)v + i [ t @ R (\\X R (s)f ml )U(t,s)F(X R (s))ds. 

Jo 

The existence and uniqueness of X R G L P (J7; C(0, T; H 1 )), the solution of 
(3.7), is easily obtained by a fixed point argument since the nonlinear term 
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we may then define a local solution X to (1.12) on a random interval 
[0, t*(v)), where t*(v) = lim^ +00 t r almost surely, by setting X(t) = X R (t) 
on [0,t r ]. Then for any stopping time r < r* we have constructed a unique 
local solution with paths a.s. in C([0, r], H 1 ). It follows from the construction 
of the stopping time r* that if r* < +oo, then limsup t _^ T , ||X(t)||jji = +oo. 
Let us now prove that if v G H 2 , then the maximal stopping time satisfies 
the following alternative: 



We note that the random propagator commutes with derivation. Hence, if 
v G H 2 , then U(-,0)v G C([0,T],M 2 ). We easily deduce, using (3.1) and in- 
terpolating M 1 between H 2 and L 2 , that U(-,0)v G C^([0, T],!! 1 ) for p G 
[0,1/4). By a fixed point argument in H 2 and equation (3.7), we conclude 
that X G C^([0,r],lH 1 ) for any stopping time r < t* and for the same max- 
imal time existence r*. Hence, using the condition on r* and uniform con- 
tinuity of X in H 1 , we get that (3.8) holds. 

Remark 3.1. We were not able to prove the global well-posedness for 
(1.12). Due to the lack of Strichartz estimates, we cannot control the evo- 
lution of the H 1 norm. Even though the deterministic energy provides a 
control on the H 1 norm because we are in the subcritical case, its evolution 
for a solution of (1.12), which is given in the next lemma, involves terms 
which are not well controlled. However, we cannot really conclude to the 
real occurrence of blow up or not in this model. It is clear that on a physical 
point of view such a phenomenon should not occur. 

Lemma 3.1. Let the functional H be defined for u G H 1 (R) by 



(3.8) 



r* = +oo or lim ||X(t)|| H i = +oo. 




Then for any stopping time r such that r < r* , we have 
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dX 2 dX 



dx dx 



2 

X2 dx ds 



d x \X\\ 2 d x \X2\ 2 dxds. 



Proof. The first equality follows by Stratonovich differential calculus 
applied to the functional H and because the process X is the solution of 
(1.12). The calculation can be made rigorous by localization (H is C 2 but 
not bounded) and regularization through convolution. The second equality 
is obtained writing the evolution of H in its Ito formulation, that is, 



where we used the unitary of the Pauli matrices and o\. = at, for k = 1,2,3. 
Easy calculations lead to the expression given above. □ 

4. Diffusion limit of the Manakov PMD equation: Proof of Theorem 1.3. 

The aim of this part is the proof of the convergence result given in Theo- 
rem 1.3. For this purpose we have to cutoff equation (1.11) in order to get 
uniform bounds, with respect to e, of high order moments of the H 2 norm 
of the solution. Let us denote by Xf" the solution of the cutoff equation 



will enable us to prove tightness on /C. 

(2) We use the perturbed test function method to get convergence of the 
generators in some sense [17, 20, 24]. This method formally gives a candidate 
for the limit process. 

(3) Setting Zf = (X& , \\Xf- (-)Hhi), we then prove that the family of laws 
C{Z*) = Po is tight on K and we deduce that the process Z^ con- 
verges in law, up to a subsequence. 




3 
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(4) Combining the previous steps and using the martingale problem for- 
mulation, we identify the limit and conclude to the weak convergence of the 
whole sequence X £ . 

(5) Finally, we get rid of the cutoff and we conclude that the sequence 
(X £ ) £> q converges in law to X in ^(IH 1 ) using the Skorokhod theorem. 

4.1. Uniform bounds on X^. Recall that a unique solution eC(R+, 
H 3 ) of the following equation exists (see Section 2): 

(4.2) -^o-i*fd^i(t)-^o- 2 *f^ 2 (t) 

+ ©ii(ll*f(*)l&)/(«f(*))* = o. 

A solution X^ to (4.1) is then easily deduced from X^(t) = Z" 1 ^)^ (t). 

Lemma 4.1. Let ^ G H 3 and Vl/f 6e £/ie solution of (4-2); then for all 
T>0 there exists a positive constant C(R,T) independent of e, such that, 
a.s. for every t in [0,T], 

||*?(t)|| H 3 <C(R,T). 

Similar bounds hold for X^(t) = Z" 1 (t)fy^(t) for any t £ [0,T] since Z~ l is 
almost surely bounded. 

Proof. The bounds on the H 3 norm are obtained using an energy 
method. Using a regularization procedure, the Ito formula applied to 
\\d x ^f (i)||| 2 and equation (4.2), we obtain for all t £ [0,T] 

\\9x*?(t)\\h = \\dMh+2 f\d x *?(s),dd x *?(s)) 

Jo 

+ ^ fwd^^wl ds, 

£ ~ Jo 

hence, 

\\d^f(t)\\b<PxMb 

+ 2 refldi^WH^jii^/^wjii^ii^fwiiycfa 

JO 

< + C{R) f \\d x ^{s)\\l 2 ds. 

Jo 
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By the Gronwall lemma we deduce that 

\\d x *?(t)\\& < \\d x Ml^MC{R)T). 

Using the same procedure for ||(9^X^||^ 2 , the Gagliardo-Nirenberg and Young 
inequalities, 

\\dl^(t)\\l 2 -\\dlM\h 

<c/ t G fl (||^( s )||^)((ll^(«)ll^ + l)||^^( s )||2 2 
Jo 

+ \\^(s)\\U\d x ^(s)\\t 2 )ds. 

By Sobolev embeddings, properties of the cutoff function and again the 
Gronwall lemma, we conclude 

\\d^(t)\\l 2 <\\dlM\l*C(R,T). 

A bound on ||c^X^||? 2 may be obtained similarly using the previous esti- 
mates and the Gronwall lemma. □ 

Remark 4.1. To prove the convergence result, we need initial data in 
H 3 (R). We will explain later where exactly we need this extra regularity, 
but this is mainly due to the fact that we prove tightness in C([0, T], H 1 ). 

Remark 4.2. Note that we first prove convergence in law for the couple 
of random variables (X^, This is due to the fact that the cutoff 

is not continuous for the weak topology in H 1 or for the strong topology in 
Hj oc . These arguments have already been used in [11]. 

4.2. The perturbed test function method. Note that the process X^ is 
not Markov due to the presence of u £ . However, (X^, v £ ) is Markov, by con- 
struction of v. We denote by its infinitesimal generator. Let us compute 
S£ E R j for / sufficiently smooth such that / maps H _1 x § 3 into R and is of 
class C%. Let (•,•) be the duality product between H 1 and H _1 . Then, for 
e > and for X^, the solution of the Manakov PMD equation (4.1), 

f(X?(t),u e (t))-f(v>v) 

= f(X?(t), u e (t)) - f(v, v £ (t)) + f(v, v £ (t)) - f(v, y) 

= (D v f(v,v&)),X?{t)-v) + R{X?{t),v) 

+ f(v,u £ (t))-f(v,y), 

where 

R(X*(t), v) = f\l - 9)(D 2 J(v + 0{X?(t) - v))(X?{t) - v),X?(t) - v) d6 
Jo 
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and Dtf{v) £ C{M~\W). Thus 



-E(/(Xf (t),u e (t)) - f(v,y)\(X(0)M0)) = (v,y)) 



E[(D v f(v,iy e (t)), 



+ 



Xf(t)-v 



E ( R(X*(t),v) 



t 



X(0)=v) +E 



(X(0),u(p)) = (v,v) 
f(v,u £ (t))-f(v,y) 



t 



i/(0) = y). 



We know by Theorem 1.1 that if v € H 3 , then X* G C 1 ([0,T],H 1 ). Thus, by 
the mean value theorem, equation (4.1), the almost sure boundedness of u, 
Lemma 4.1 and the conservation of the L 2 norm, 

1 



t 



X?(t)-v\\ h 2 



< sup \\d 3 X?(s)\\ L 2 
se[o,t] 



< sup 

s€[0,t] 



j*(v £ (t))d x X*(s) 



+ 



< 



V do 
e 2 



+ ||e i? (||^( s )|| 2 I1 )^ (S )(^( s ))|| L2 

C(R,T) + 2RC\\v\\ h 2. 



Thus, by the boundedness of D 2 f, the continuity of t >-)■ X^(t) in L 2 and 
the previous bounds, we conclude that 

R{X % {t),V) <C(R,T,e) sup \\D 2 J( W )\\ m - 1>ml) (l + \\v\\^)\\XS)-v\\^ 

and the right-hand side above tends to zero as t goes to zero. Now, we 
perform the change of variables t' = t/e 2 to get 

±E(f(v, u e (t)) - f{v, y)H0) =y) = »(?)) - f(v, y)\u(0) = y). 

Thus, using the Markov property of the process v, and using (4.1) again, we 
get an expression of the infinitesimal generator of the Markov process 



J? E H f(v, y) = Km -(E(f(X?(t), u £ {t)) - f(v, y)\(X(0), u(0)) = (v, y))) 



{D v f(v, y),d t X^t)\ t=0 ) + -=Jf v f(y, y) 
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(4-3) 

= (D v f(v,y), ^ + &r(\M&Fv («)) 

1 / (9?; \ 1 

- -(D v f{v,y),b'*{y) —J + ^JSf„/(t;,y), 

where jSf^ is the infinitesimal generator of v and <&l v its domain. The per- 
turbed test function method gives (by identifying its infinitesimal generator) 
an idea of the limit law of the sequence (X r ) £> q. It provides in addition con- 
vergences that are useful to prove the weak convergence of the sequence of 
measures (C(X r )) £> q. 



Proposition 4.1 (Perturbed test function method). There exists a lim- 
iting infinitesimal generator (J? R ,S> R ) such that for all sufficiently smooth 
and real-valued functions f € 3> R and for all positive e, there exists a test 
function f £ and positive constants C±(K) and Ci(K) satisfying 

(4.4) sup \f £ (v,y)-f(v)\<eC 1 (K), 

veB(K) 

yes 3 

(4.5) sup \J? e R f £ (v,y)-J? R f(v)\<eC 2 (K), 

veB(K) 

where B(K) denotes the closed ball o/H 3 (M) with radius K. 



Proof. The idea is to prove that for all suitable test functions /, one 
can find a function f £ of the form 

(4.6) f £ (v, y) = f(v) + efiv, y) + s 2 f(v, y), 

such that Proposition 4.1 holds. We plug this expression of f £ into (4.3) and 
formally compute the expression of ££ R f £ : 

^ £ R fe(v,y) = (^D v f(v), l A^ + l e R (\\v\\ 2 ml )F y (v) 
D v f\v,y),b'cr(y)^- 
+ ^f 2 (v,y) + -2> v f\v,y) - -/D v f(v),b'*(y) — 

(4-7) 

+ e/D v f\v, y), ^ + i& R (\\vf ml )F y (v) 
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-e^D v f 2 (v,y),b'cT(y)^ 

+ e 2 (D v f\v, y), ^° + ie fl (||t;|| ^)JFJ, («) 

and we notice that ££ v j{y) is identically zero because / does not depend on 
v = iyx^v-i). The aim is to wisely choose the functions f 1 and f 2 and the 
regularity of / so that -S^/e is wen defined and that f e and S£^j e converge 
in the sense of Proposition 4.1. In particular, we need to cancel the terms 
with a factor 1/e and we need the terms with factors e or e 2 to be 0(e) on 
bounded sets. In order to cancel the 1/e terms, we look for a function f 1 
solution of the Poisson equation 

(4.8) J?„f\v,y) = /D v f(v),b'v(y)^ 

By Corollary 5.1, we know that 

E A ( 5 ») = Vj = l,2,3. 

We deduce that (D v f{v),b'a(y) ff), which is a linear combination of m, = 
gj{y) [see (1.4)], is of null mass with respect to the invariant measure A. 
Hence, (D v f(v),b'a(y) is a function of y G S 3 , which satisfies the as- 
sumptions of Proposition 5.1, provided that / is sufficiently smooth, that 
is, / G C 1 (IH _1 ) and dSL 2 . It follows that the solution f 1 of the Poisson 
equation (4.8) can be written as 



(4.9) 



where 



f\v, y) = J?~ l (I D v f(v), b'a(-) 



(y) 



D v f(v),b'a(y)^- 



r+oo 

(4.10) *(v) = E(-*(y(t))\u(0)=y)dt. 

Jo 

By Proposition 5.1, there is a positive constant M such that 

(4.11) ||£(i/)||oo<Af VyGS 3 , 

and f l {v,y) is a continuous bounded function of y for i> G L 2 . We now have 
to choose the function f 2 , but we cannot choose J£ u f 2 cancelling the terms 
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because they do not satisfy the null mass condition with respect to A. Hence, 
we look for a solution f 2 of the Poisson equation 

^uf(v,y) = -(D v f(v),iQ R (\\v\\ 2 m )F y (v)) 
+ (D v f(v),iQ R (\\v\\ 2 m )F(v)) 

(412) ' dv 

D v f(v,y),b'<r(y) — 

where, due to (4.9), 

D v f\v,y),b'a(y) — 



(4-13) =-{b'Y{Dif(v)a{y)-My)Q- x 

Moreover, thanks to expression (4.13), the Fubini theorem and Corollary 5.1, 
-m A (^\v,y),b'v(y)^\\ 



(4.14) 



A(9k(Ht))9M0)))dt 



E A (g k (u(t)) 9j (u(0)))dt 



fc=i 



where 7 = (b') 2 /6~/ c . Provided that / is of class C 2 (M.~ 1 ) and wei 1 and 
because f l (v,-) is of class C|(S 3 ) for any u € H 1 , we can now define, by 
Proposition 5.1, a unique solution, up to a constant, to the Poisson equation 
(4.12). This solution f 2 is expressed as 

f(v,y) =J?- 1 ((D v f(v),i@ R (\\v\\ 2 ml )(F y (v) -F(v)))) 
-^(/ Dv f\v,y),b'cT( y )^ 
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(4.15) 



(D v f(v),ie R (\\v\\ 2 ml )F(v,y)) 
3 



k,l=l 



where 



D v f{vUb'fa{y)^, 



n+OO 

F(v,y)= I E(F u{t) (v) - F(v) \ u(Q) = y) dt 



and 

f fcl i(y) = y o ( y E(»fc(f (o) = y)ds - «ft 



+oo / f+OO ^ 



and 

»+oo / r+oo 



a{y) = J (/ E( < r( I /( S )) < r( I /(t))|f (0) = y) ds - ^ J dt. 

Replacing Jjf„f and J£ v f in (4.7), respectively, by the right-hand side of 
(4.8) and (4.12), and using expression (4.14), we get 

*e*fe(v,V) = (D v f(v), + y) + ^(|b||^)^) 

fe=i N 



(4-16) + e (^/ 1 K 2 /),^0 + iG /? (lK 1 )^K 



and we define the limiting operator by 



(4.17) 



7 V"^ / n 2,/ \ 9w 

fc=l x 
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Hence, if we define & R as the space of functions which are the restriction 
to H 3 of functions / from H" 1 into R of class C 3 (H~ 1 ) and such that / 
and its first three derivatives are bounded on bounded sets of HI -1 , then the 
functions f l and f 2 are well defined for / E 3l R . Moreover, if / E @ R , then 
S£ R f e is well defined for v E H 3 . 
We now write that 

sup \f e (v,y)- f(v)\<e sup \f x (v,y)\+e 2 sup \f 2 (v,y)\ 

v£B(K) v&B(K) v&B(K) 



y& 



yes 3 



yes 3 



and use the following result, which is proved in Section 6. 

Lemma 4.2. Let f E 3> R and f 1 and f 2 be, respectively, solutions of 
(4.8) and (4.12). Then 

sup l/V^CiW and sup \f\v,y)\ <C 2 {K). 

veB(K) veB(K) 

ye§ 3 yes 3 

This proves the first convergence of Proposition 4.1. With 5£ R f(v) given 
by (4.17), the second convergence (4.5) in Proposition 4.1 follows from (4.16) 
and the next lemma, which is proved in Section 6. □ 

Lemma 4.3. Let f E $> R and f 1 , f 2 be, respectively, solutions of (4-8) 
and (4-12). Then 



sup 

veB(K) 
yes 3 



D v fHv,y), ^ +i@R (\\v\\ 2 ml )F y (v) 



sup 

v£B(K) 
j/<E§ 3 



D v f 2 (v,y),b'a(y) — 



sup 

v£B(K) 
3/GS 3 



D v f(v,y)^^+i& R (\\v\\ 2 nl )F y (v] 



<Ci(K), 
<C 2 (K), 
< C 3 (K). 



4.3. Tightness of the family of probability measures (C(Z r )) £> q. To prove 
tightness on K, of the sequence of probability measure C(Z R ) = Po (Z R )~ l , 
we need to obtain uniform bounds in e on Z R in the space 

(C([0,T],H 2 )nC a ([0,T],I- 1 )) x C 5 ([0,T],R) 

for suitable a, 5 > 0. Note that uniform bounds of X R in C([0, T], H 2 ) are 
given by Lemma 4.1. The perturbed test function method will enable us 
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to get the uniform bound in C Q ([0,T],H -1 ). Such bounds cannot be di- 
rectly obtained using (4.1) because of the 1/e term. In order to obtain such 
bounds, we use again the perturbed test function method for convenient test 
functions. Let (ej)j^* be a complete orthonormal system in L 2 . Recall that 
(•, •) is the duality product between EI 1 -!! -1 and (•, -)j2 the inner product in 
L 2 . By definition of H s ,sSK, we can define a complete orthonormal system 
(ej)jgN* on H 1 from (ej)jen*> 



i«iii-i = ii(i+e a r 1/a ^ 



L 2 



where ej = J r-1 ((l + £ 2 ) -1 / 2 £f) for any j £ N*. We denote by (fj)jeN* the 
family of test functions in Qi R defined by 

Vh ^ fj( v ) = (ej,v). 

For w£H 3 , we also consider particular perturbed test functions fj >E of the 
form 

(4.18) f jje (v,y) = f j (v)+efj(v,y), 

where, for all j in N*, fj(v,y) = (ej,ip (v,y)) for a given function tp 1 with 
values in H 2 . We now choose if 1 as a solution of the Poisson equation in y: 

(4.19) ^ l p 1 (v,y)-b'a(y) — =0, 
whose explicit formulation is given by (see Proposition 5.1) 

(4.20) *\v,y) = -}f<T{y)^, 

where <x(y) is given by (4.10). We point out that ip behaves in its first 
variable like and is linear in v. Consequently, for all j in N*, 



dx 



ei) ^^^ + ^(||Xf(i)|| 2 1 )^ e(t) (Xf(t)) 
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/ 9 fido d 2 X R (t) 

+ ie R (\\x*(t)fwi)Fu e (t)(x?(t)) 

For all t G [0, T], we define the process M R with values in HI" 1 given for any 
j in N* by 

= f jj£ (X R ,u £ )(t)-f j>£ (v,y)- [ J? e R f h£ (X*{s)Ms))ds 



(e j> X R -v)+e(e j ,i P 1 (X R ,u £ ) - ^{y,y)) 

- f ^ £ R f h£ {X^)Ms))ds. 
Jo 



Given the fact that ^£ R is the infinitesimal generator of the continuous 
Markov process (X R ,u £ ) and Jzf/Vj.e is well defined because fj G @ R , then 
(ej, M R (t)) is a real- valued continuous martingale. Moreover, it is a square 
integrable martingale, as follows from the bounds on the H 3 norm of X R 
obtained in Lemma 4.1. To prove tightness of the family of probability mea- 
sures C{Z R ) on /C, we need estimates of moments on the processes X R and 
ll^e^Olle 1 ' Before proving these estimates we introduce a process Y £ R close 
in probability to X R for which it will be easier to get those estimates, us- 
ing, in particular, the Kolmogorov criterion. The idea is to use Lemma 4.4 
below to get tightness of the family C(Z R ) from convergence in law of a 
subsequence of Y R . 

Lemma 4.4. Let us define the process Y R as 
(4.22) X R (t)-Y R (t)=e^ 1 (v,y)-^(X R (t),u £ (t))) VtG[0,T]; 
then for all 5 > 0, 

P(\\X R - Y R \\ cmT]m >5)< |d(r, R). 

Proof. Using the Markov inequality and Lemma 4.1, we get for all 

5>0, 

P( sup \\X R (t)-Y R (t)\\ ml >5) 

V t£[0,T] 7 

<^E( sup \\ip\X R (t),v £ {t))-V X {v,y)\\ m i 
d v ie[o,T] 
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< |E I sup 



5 



te[o,T] 



<£™C(T,fl), 
where M is given by (4.11). □ 

Note that the process Y R is also defined by the identity, for all j in N*, 
fo, = (ej,X? (t)> - e(e,V(«, V) " "«(*))> 

(4.23) = < ei ,M«(t)) + (e^) + f J? £ R f h£ (X?(s)Ms))ds 

Jo 

yt £ [0, T]. 

Lemma 4.5. For aZZ 1 > e > 0, i/iere exist three positive constants C\(T, 
R), C2(T,R) and C%{T,R) depending on final time T and on the cutoff 
radius R, but independent of e, such that 

(4-24) ^\YX\fc (mm2) ) KCifrR), 

(4-25) Edll^llca^ri^-!)) <C 2 (T,i2), 

(4-26) E(||||y/||^|| c , ([0)T]iR) )<C 3 (T,i?), 

where < a < \ and 5 = a/3 > 0. 

Proof. Thanks to Lemma 4.1, we know that the solution X R of (4.1) is 
uniformly bounded, for all e, in H 3 by a constant C depending on R and T. 
We conclude, using the explicit formulation of ip\ given by (4.20) and (4.22), 
that (4.24) holds. 

To prove inequality (4.25), we first need an intermediate estimate that 
will be proved in Section 6. 

Lemma 4.6. There exists a positive constant C(R,T) such that for all 

t,se[0,T] 

EQlY^-Y^f^KC^TXt-s) 2 . 
Then we deduce from Lemma 4.6 

E(||'5^ R |Iwt.4([o,T],H- 1 )) — 

for any 7 < 1/2. We use the Sobolev embedding W 7,4 ( [0, T] , H^ 1 ) <-> C a ([0,T], 
M^ 1 ) for 7 — a > 1/4 and 7 < 1/2, which implies a < 1/4. Thus, we deduce 
the second inequality (4.25). 
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It remains to prove the last bound (4.26). Note that for t,s € [0, T] 

\\\Y e R (t)\\^-\\Y e R (s)\\U 

<C sup \\Y R (r)\\ m \\Y R (t) -Y R (s)\\ ml 

re[0,T] 

<C sup \\Y R (r)\\ m sup Hl^Wdl^W-^Wtfi- 

re[0,T] r£[0,T] 



It follows that if 5 = a/3, 

^ n D ,m ilv-R^Ml 5 / 3 [n 

l e Nc^GO.T],]!- 1 )" 



m fl (oiiHiiic«(p I Ti«<c sup iin a (r)iiSfnn Ji||1/8 



re[0,T] 

Inequality (4.26) is then implied by the Holder inequality, (4.24) and (4.25). 

□ 

Remark 4.3. The extra H 3 regularity is needed precisely in the first 
step of the above proof in order to estimate the H 2 norm of Y R , which 
involves the gradient of X R . 

Proposition 4.2. The family of laws (C(Z r )) £> q is tight on JC. 

PROOF. We set Z R = (Y R , HF/OHhi)- Denoting by B{K) the closed 

ball of (C([0,T];H 2 (M))nC Q ([0,T];H~ 1 (^))) x C 5 ([0, T]; M) with radius K, 
for a and 5 as in Lemma 4.5, we deduce using the Ascoli-Arzela and Banach- 
Alaoglu theorems that B(K) is compact in /C. Using the Markov inequality 
and Lemma 4.5, we get 

¥{Z R iB{K)) 

- ^^("^^{ll^llcdO.T]^ 2 ), 1 1 1 1 C« ([0,^] ,-H- 1 ) ; IHI^IIlIillc^ICT])}) 

< 1 max(C 1 1/4 (T, R) , C 2 (T, R) , C 3 (T, R) ) . 

We conclude that the family of laws (C(Z r )) £> q is tight on K, and by the 
Prokhorov theorem we obtain the relative compactness of the sequence of 
laws {C(Z R )) £>0 , that is, up to a subsequence, the sequence C{Z R ) weakly 
converges to a probability measure C(Z R ) where Z R = (X R ,j R ). We may 
now use Lemma 4.4 to prove that the family of laws C(Z R ) is tight. Indeed, 
it easily follows from Lemma 4.4 and the above convergence in law that for 
all g € C b {K) 

KmE(g(Z?))=E(g(Z R )). n 
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4.4. Convergence in law of the process X R . In order to get the con- 
vergence in law of the whole sequence (X r ) £> q, it remains to characterize 
the limit, that is, to prove that X R = X R , the solution of (3.6), and that 
j R (t) = HX^t)^ for any t G [0,T]. The tool here will be the use of the 
martingale problem formulation introduced by Stroock and Varadhan in [26]. 

Proposition 4.3. The whole sequence X R converges in law to X R in 
C{%T]M). 

Proof. In order to prove that any subsequence of X R converges to 
the same limit X R , the solution of (3.6), we will prove the convergence of 
the martingale problem for suitable test functions / G @ R . To this purpose 
let us define, for a G H 1 with compact support, the particular test function 
f a (-) = (a,-), so that f a G 3> R . From this particular choice, we construct a 
perturbed test function / 0)£ , 

fa,e( v ^) = fa(v) + ef^(v,y) + e 2 f%(v,y), 

obtained thanks to Proposition 4.1. The correctors /* and f% are chosen 
to be the solution of the Poisson equations (4.8) and (4.12) for f a . Let us 
denote by Z R a subsequence converging to Z R and define the EI -1 valued 
process Nf (Z R (t)), associated to (4.1), 

(a,N R (Z R (t))} = f a , £ (X R (t),u £ (t)) - f a>e (v,y) 

- [ ^ R f a!e (X R (s),u E (s))ds, 
Jo 

where S£ R is given by (4.7). We also define the process N R (Z R (t)), 

(a,N R (Z R (t))) = f a {X R {t)) - f a {v) - f J? R f a (X R (s))ds, 

Jo 

where 5£ R is given by expression (4.17). Moreover, we denote by ^£ R R the 
operator whose expression is given by (4.17) replacing ||X R (£)|| H i by r y R (t) 
in the cutoff function. Let us now define (a,N R (Z R (t))} by 

(4.27) (a,N R (Z R (t))) = f a (X R (t))-f a (v)- f J? R f a (X R (s)) ds. 

Jo 

The process (a, N R (Z R (t))) is a real continuous martingale because (X R , u £ ) 
is a Markov process and because £^ R f a ^ £ is well defined since X R (t) el 3 . 
Moreover, it is a square integrable martingale, as follows from the bounds on 
the H 3 norm of X R obtained in Lemma 4.1. The above martingale property 
implies that for all t, s G [0, T], t > s, 

E[(a,N R (Z R (t))-N R (Z R (s)))\a(Z R (u),u £ (u)),u<s] = 0. 
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It follows, in particular, that for all test functions hi,..., h m G C&(Hj oc x R) 
and < ti < ■ ■ ■ < t m < s < t, 

m 

E (a,N«(Zf(i))-Nf(^( S )))n^(^(^)) =0- 

Using Proposition 4.1, Lemma 4.1 and the boundedness of the functions /ij, 
we get 

/ m n 

E (a,N*(Z*(t)) - N R (Z R (t)) - Nf + N«(Zf (,))) n^(^(^)) 

V 3=1 / 

<eC(R,T). 

Let us consider a cutoff function x_R £ C^°(^) satisfying 

1, ]fueBjc(Ro), 
0, ifn^^(2ii ), 



Xi?o(«) 



where Bic(Ro) denotes the closed ball of radius i?o of the space K, and 
Rq is chosen such that X R G B/c(Ro) a.s. (see Lemma 4.1). Note that by 
continuity of the functions XRo an d {^j}je{i,...,m}; respectively, in fC and 
x R, by continuity of / o (0 for the weak topology in H 1 , by continuity 
and boundedness of @r in C([0,T];R), by continuity of F from H 1 to H _1 
and the bounds on F(X R (t)) obtained thanks to Lemma 4.1, the function 

m 

(a^Z^mXRoiZ^Uh^it,)) 

3=1 

is a bounded and continuous function of Z R from /C into R. We deduce by 
convergence m law of Z R to Z R in K, since the test function a is compactly 
supported, that for all t, s G [0, T], t > s 

(4.28) E^a,^^^)-^^^)))^^^)!!^^^)^ =0. 

Since, almost surely, X R belongs to the closed ball B/c{Rq), we deduce 
that almost surely X R G B)c(Rq). Thus, we conclude from (4.28) that (a, 
N R (Z R (-))) is a continuous square integrable martingale with respect to the 
filtration Qt = a(Z R (s), s < t) and this holds for any a G H 1 with compact 
support. 

In order to identify the equation satisfied by X R , we consider, for a, b G 
H 1 with compact support, the function g a ,b(v) = fa( v )fb(v) G ^ R and the 
perturbed test function g a ,b,e, 

9a,b,e{ v ^y)=9a,b{v)+eg l ah (v,y)+e 2 g 2 ah (v,y), 
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obtained thanks to Proposition 4.1. Thus, functions g\b{v,y) and g^ b (v,y) 
are chosen to be solutions of the Poisson equations (4.8) and (4.12) for g ab . 
Let us now define the real-valued continuous martingale 

H^ b , £ (Z R (t))=g aA£ (X R (t),u £ (t))-g aA£ (v,y) 

- f J? R g a , b , £ (X R ( S ),v £ (s))d S . 
Jo 

Using the same arguments as before, we may prove that 

limE^(H^ e (Zf(t)) - Hj^^W^XBoC^n^C^i))^ 

= E ^(H. R b (Z R (t)) — Ji^ b (Z R (s)))xR (Z R )Y\^hj(Z R (tj)^ , 

where 

n R b (Z(t)) = g a , b (X R (t))-g a , b (v) 

1 'j? R R g aih {X R {s))ds. 

From the above convergence and the martingale property of H R b e (Z R (t)), 

we deduce that H R b (Z R (-)) is a continuous real-valued martingale. A clas- 
sical computation then shows that the quadratic variation of the martingale 
N R (Z R {t)) defined in (4.27) is given by 

(b,((N R (Z R (t))))a) 

rt J? R (f a (Z R (s))f b (Z R (s))) - f a (Z R (s))J? R f b (Z R (s)) 



-MZ R (s))J? R UZ R (s))ds. 

Applying the operator ££ R R , respectively, to the test functions f a and g a ^ b , 
we obtain that 



Sf$fa(Z"(t)) = ( a, I + -± \ + ie R (\\^(t)r m )F(X J 



and 



^i«a, t (Z"(f)) 
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fc=i x ' x 

We deduce that the quadratic variation is given by formula (3.5) with X 
replaced by X R . Thus, using the martingale representation theorem, we can 
write the ^-martingale N R (Z R (t)) as the stochastic integral 



(a,N R (Z R (t))) = ^y /*j7a,<7* 



where W= (Wi,W2,W3) is a real- valued Brownian motion on a possibly 
enlarged space (Q,Q,Q t ,F). We deduce that (.X^, W) is a weak solution in 
C([0,T];Ul c ) n CuQO.TIjM 1 ) nL^(0,T;EI 2 ) of the equation 



(4.29) 



idX R {t) + (^%X R (t) + e R ( 7 *(i))F(X*(t)) ) cit 

3 

+ iy/T%2<T k d x X R (t) O dW k (t) = 0, 



fc=l 



I = "£i 3 . 



The next step consists in proving that almost surely j R (t) = \\X R (t)\\^ s i. 
Using the Skorokhod representation theorem, we can construct new random 
variables (that we still denote Z R , Z R ) on a new common probability space 
(fi, J 7 , Tu IP) with, respectively, C(Z R ) and C(Z R ) as probability measures 
and with values in K, such that 

\miZ R = Z R , P-a.s. in /C. 

Since X R G L°°(0, T; M 2 ), we deduce using (4.29) that X R G C([0,T];L 2 ). 
Hence, applying the Ito formula, it is easy to see, since Qr is a real-valued 
function, that almost surely 

\\X R (t)\\ L2 = 1Mb = \\X R (t)\\ L 2 Vt G [0,T],We > 0. 

Thus, we deduce the strong convergence of X R (t) to X R (t) in L 2 , a.s. for 

each te [0,T]. Since X R converges to X R in L~(0, T; H 2 ), we get using 
Lemma 4.1 that 

l|£ R |U-(0,T;H3) < Hminf H^Tf |U« (0 ,t;BP) < C(R,T), P-a.s. 
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Interpolating H 1 between L 2 and H 2 , we conclude that 

(4.30) \im\\X R (t)-X R (t)\\ m i=0 Vte[0,T], P-a.s., 

and X R G CQO.T];!! 1 ); it follows that, almost surely for all t in [0,T], 
-f R (t) = \\X R (t)\\l ll and X R is a solution of (3.6). Thus, the limit in law 
of X R is unique and is given by the solution X R of (3.6). 

The final step consists in recovering the convergence in law in C{ [0, T] , H 1 ) . 
Since Y R is uniformly bounded in e in C a ([0,T],M~ 1 ) n C([0,T];H 2 ) with 
< a < 1/2, we deduce that it is a.s. uniformly bounded in e in 
with (3 = a/3. Moreover, using pointwise convergence (4.30), expression 
(4.22) and uniform bounds (4.1), we get pointwise convergence in H 1 of 
Y R to X R . We conclude that Y R converges in law to X R in C([0,T],1[ 1 (M)) 
and by Lemma 4.4, the convergence in law of X R to X R in C([0, T], H 1 (R)) 
follows. □ 

Remark 4.4. Using the Arzela-Ascoli and Banach-Alaoglu theorems, 
Lemma 4.5 and the Tychonov theorem, we deduce that (C(X r ))ji^ is tight 
on /C N . Thus, the same arguments as above lead to the convergence in law 
of (X R ) Rm to (X R ) Ren (see [11]). 

4.5. Convergence of (X e ) e> Q to X. Using the Skorokhod theorem, we 
can construct new random variables X R , X R on a common probability space 
(tt,T,y t ,F) and with values in C([0, T], M 1 ) such that for any R > 0, 

f n R = n R ~ 

\% Me D ' and X? -»• X R , P-a.s. in CdO^},!! 1 ). 
[fl R = H R , e^o 

We define the escape times t r and t r associated to the cutoff: 

r R = inf{te[0,T},\\X R (t)\\ ml >R} 

and 

r £ R = M{te{0,T],\\X R (t)\\ M i > R}. 

Let X £ and X be the processes, with values in £'(E[ 1 ), defined, respectively, 
by X £ (t) = X R {t) for t < t r and X(t) = X R (t) for t < t r , X(t) = A for 
t > t* = hniR_j. +00 t r . Then if r < r* a.s. is a stopping time, the process 
X £ converges to X a.s. in C([0, r], H 1 (M)). Hence, the convergence in law in 
^(IH 1 ) follows. 

5. Study of the driving process v. We recall in this appendix some re- 
sults obtained in [18, 22] about the driving process v. 

Proposition 5.1. The process v = (^1,^2)* is a Feller process that evolves 
on the unit sphere S 3 of C 2 ~ M 4 . Furthermore, it admits a unique invariant 
measure A, which is the uniform measure on § 3 , under which it is ergodic. 
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For all f € C 2 (§> 3 ) satisfying the Fredholm alternative (or null mass condi- 
tion) K\(f(u)) = Jg 3 f(y)A(dy) = 0, the Poisson equation J£ v u(y) + f(y) = 
admits a unique solution of class C 2 (§ 3 ), up to a constant, which can be 

written as u(y) = J +oc E[/(i/(t))|i/o = y] dt. 

Let us recall that <r(v(t)) = a\m\ J rGi'm<i + 0-31713 where rrij{t) = gj{v{t)). 
We now state a result related to the effect of the random PMD on the pulse 
evolution. 

Corollary 5.1. (1) The process m = (mi, 777.2, 7713) € § 3 is a Feller pro- 
cess with a unique invariant measure Ao g^ 1 under which it is ergodic. 

(2) For j = 1,2,3, E A (s») 
Sjk/3. As a consequence, 



E AoM™) =0 andE A ( gj (u(t))g k {u(t))) 



E A (iV(X)) = |(2|X 2 | 2 - {X^Xi, 



3^1^21 

E A (N 2i „(X)) = l(2\X 1 \ 2 
(3) For j, k = 1,2, 3, 



|Xa I )Xo. 



f 

Jo 



E A [ 9j (u(0))g k (u(t))]dt 




ifj = k, 
ifj k, 



where 7 C is the constant appearing in (1.8). 
6. Proof of technical lemmas. 



PROOF of Lemma 4.2. Let v be in H 3 . Using the explicit representation 
(4.9) of f 1 , we obtain, since D v f(v) G H 1 (R), that 



D v f(v),b'a(y)^- 



<b'\\D v f(v)\\ m i 



dv 



dx 



+00 



E(g j (v(t)M0) = y)dt 



Moreover, by Proposition 5.1 the integral J +o °E(g , j(^(t))|z/(0) = y)dt con- 
verges because gj is a bounded function of v € S 3 . Since u h-» D v f(v) is a 
continuous function which is bounded on bounded sets of EI -1 , we deduce 
that 

sup \fHv,y)\<b'C(K). 

veB(K) 

2/es 3 

The function f 2 given by (4.15) may be bounded using the same arguments. 
Indeed, 

(D v f(v),iQ R (\\v\\ 2 m )F(v,y)) < \\D v f(v) 



\F(v, 
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Since for all v G H 3 , y i— > F y (v) — F(v) is a function of class C 2 on § 3 , with 
values^in HI -1 , satisfying the null mass condition of Proposition 5.1, the 
term F(v,y) is bounded. Moreover, v i— > F y (y) — F(v) is bounded in 
on bounded sets of H 1 by the continuous embeddings H 1 (R) L 4 (M) and 
L 4 / 3 fR) ^->M -1 flRV In addition, 



kl=l 



< 



k,l=l 



n2 ,, <9v <9u 



+ 



D v f(v),a k ai 



8 2 v \ 
dx 2 1 



< 



C ||^/W|U (H - llH i 



<9r 



<9.r 



+ I|A,/Mlk 



d 2 v 



dx 2 



Since v h-> D v f(v) and t> i— >• D 2 f(v) are bounded on bounded sets of EI 
we conclude the proof of the lemma. □ 



Proof of Lemma 4.3. Replacing f 1 by its expression (4.9), we get 
'D v f\v,y), 



idn d 2 v „ ... „« . „ . . 



Dlf{v)b'a(y) — — + l Q R {\\v\\l 1 )F y {v) 



By the assumptions on /, v i— )• D v f(v) and u i— )• D 2 f(v) are continuous 
bounded functions on bounded sets of H _1 (R). Moreover, D 2 f(v) G £(H _1 , 
M 1 ), £)„/(«) G M 1 and G L 2 . Using the bound (4.11), we deduce that 



sup 

yGS 3 



D^yJ^g +^(11^11^)^^) 



<C(iT). 



Let us now compute the first derivative of f 2 using expression (4.15); for all 
h in H 1 and u in H 3 , 

(D v f(v,y),h) 

= (D 2 J(v)h,iQ R (\\v\\ 2 m )F(v,y)) 

+ (D v f(v),2iQ' R (\\v\\ 2 ml )(v, h) u iF(v, y)) 

+ (D v f(v),iQ R (\\v\\ 2 m )D v F(v,y).h) 



DIFFUSION APPROXIMATION THEOREM 41 

k,i=i ^ ^ 

AJM,(6O 2 ^)0). 

Taking, respectively, /t = f £ + *'©fl(||«|lHi)^( u ) and h = b>(T (y) we 
conclude 



sup 

veB(K) 

y&S 3 



<C(iT) 



and 



sup 

yes 3 



D v f\v,y),b'a{y) 



dv 
Or 



<C(K), 



since t> i— >■ D%f(v) is bounded on the bounded set of EI with values in 



^(H- 1 ,^) and & G HT 1 . □ 



Proof of Lemma 4.6. Let us recall that the family {e.;}j g N* denotes a 
complete orthonormal system of H 1 constructed from a complete orthonor- 
mal system {ej}j G pj* in L 2 and (•, •) is the duality product between H — H . 
Then 

C+oo 2 

||y«(t)-y«(,)||^ 1 = |^(e i ,y«( i )-y«( s )) 2 | . 

Using twice the Young inequality and the expression of given by (4.23) 
and (4.21), we obtain 

nn*(*)-n H wiiH-i 



d -lJ s dlx^t')dt' 



c 



e R {\\x?{t>)\^)F Ve{tl) {x?{t'))dt 
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(b') 2 a(u £ (t'))<r(u £ (t'))dlX^t')d/ 
+00 1 2 

i=l 



+ Gr(\\x b (i!)\&)f mv) (x*V)) )dt' 



We bound each term separately. Using Lemma 4.1, 



dx 2 



dt' 



<C(R,T)(t-s) 



Using that F is cubic and Lemma 4.1, 

*e*(||x.(OII*)JW*eV))<ft' 

and using Lemma 4.1 and the bound (4.11), 



dx 2 



dt' 



<C(R,T)(t- s) 



<C(R,T)(t-s) 



Finally, we bound the e 4 term that is well defined because Xf' has values 
in H 3 . Using the Cauchy-Schwarz inequality, Lemma 4.1 and (4.11), we get 
for all e < 1, 



< e 4 (6') 4 M 4 



f y + ^(IWI^W^V))) eft' 



<C(i?,T)(i-s) 4 . 

Taking the expectation and adding the previous estimates, we deduce that 

E(nn H (*)-n H wiiH-o 

< C(R, T)(t - S ) 4 + CE ( [ J] fo, M*(t) - M*( 5 )> 5 

In order to prove a uniform bound, with respect to e, of the second term, 
we will use the Burkholder-Davis-Gundy inequality and, consequently, we 
have to compute the quadratic variation ((M^(t))) of Mf (t) defined, for all 
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j G N*, by 



where Jz? £ R fj^ £ (X R (s),v £ (s)) is given by (4.21). The next lemma states that 
the process ((M R (t))) can be expressed only in terms of the infinitesimal 
generator of the Markov process v. 

Lemma 6.1. For all j in N* 

MM*(t)))ej) 

- 2(b') 2 J*/ ej , ^^))U^M^e{s)) ^{s)\ ds. 

Thus, using the Burkholder-Davis-Gundy inequality, 

E^^( ej ,M R (t)-M £ R (s))^ ^j<C(R,T)\t-s\ 2 , 

thanks to Lemma 4.1 and Proposition 5.1. Adding the previous estimates, 

EiWY^-Y^W^^CiR^t-sl 2 , 
and Lemma 4.6 is proved. □ 

Proof of Lemma 6.1. A classical computation shows that for all 



(e i ,((Aff(t)»e i >= / J? £ R (f j>£ (X?(s)Ms)) 2 )ds 
Jo 

Jo 

Now, for all j in N*, 
(/,- £ (X*( S ),z, £ ( S ^ 

I f)X R ^ 2 

+ (b') 2 e 2 L J ,*(Ms))^ r ( S ) 
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Thus, we get 

J? £ R (f jt£ (X?(s),u £ (s))f 



2< Ci , jrfc^^cy, ^5. + ie«(||jr 4 f (^[^^(^(jsrf (*)) 

- 2b'e(e 3 , ^ + iB R (\\X* (s)\\ 2 nl )F Ms) (X? (,)) 

x /e^5(i/ e (s))^_(a)\ 

x (e^M*)) A ^ ^|Sf) + l e R (\\X^s)\g 1 )F Ms) (X^s)) 



+ (6') 2 ^(^,^ £ ( S ))^( S )^ 2 

+ 2(60V( e ,,5(, £ ( s ))^ (s )^ ei)? (, £ ( S ))f ^) 
+ 2(6') 2 e 2 ^e J ,ar(^( s ))^( s ; 

x^ ei) z5F(^( s ))e fl (||xf( s )||^)^( s) (^(«))) 

The same kinds of computations for the term 2fj^J£ £ R fj^ £ lead to the result. 
□ 

Acknowledgments. The authors wish to thank the referees for fruitful 
comments. 
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